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We derive an explict bias dependent expression for electron transfer reaction rate 
from a solvated redox to a electrode through a bridged molecule of arbitrary length. 
The interaction of the solvated redox with the solvent is modelled as a classical har- 
monic oscillator bath. The effect of competing process, namely resonance tunneling 
between redox and bridge and the solvation of the redox is investigated. Plots were 
produced for the case of 5 atom bridge. Our analysis shows that for certain suitable 
value of nearest-neighbour coupling, it is possible to block electron transfer in certain 
intermediate voltage regime 

Abstract 



1 Introduction 



Electron transfer at molecular levels and , especially, between chemically active 
species has been an important and interesting field of active research for a couple of 
decades [T] [2] [31 H] • A typical quantity of interest in chemistry is the rate of transfer 
between two chemiactive species, that is, between donor and acceptor species which 
in general are solvated. Additionally, with the advent of nanotechnology and molecu- 
lar electronics, understanding electron transfer rate in molecular chains became a key 
theoretical interest. A general setup in molecular electron transfer is to study, both in 
theory as well as in actual experiment, the current-voltage response ontained on pass- 
ing current between two molecular mesoscopic junctions connected by a single long 
molecule or a chain of repeated molecular or atomic units. The important quantity of 
interest in the above scenario is the molecular conduction. The problem of obtaining 
the conduction of a such a system is well studied one. Formal expressions and rela- 
tionships like the landauer formula [S] El El an d it's several variants are available to 
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calculate the conductance of such systems [8] [9] [10] [TT] . Since conductance is due to 
transfer of electrons, the question of relationship between conductance and electron 
transfer rate was answered by A. Nitzan and co-workers |12l 113] . Nitzan derived a 
relationship between electron transfer rate (k) and conductance at certain regime. 



The rate expression for electron transfer has been derived by several methods, the 
earliest dating back to the superexchange method proposed initially by McConnell 
[141 1151 1161 1171 1181 1191 I20| . Superexchange treatment is a viable method when the 
transfer mechanism is tunneling dominated. Another plausible mechanism for transfer 
is sequential hopping. Typical cases of sequential hopping occurs when there are asym- 
metries or irregularties in the bridge connecting the acceptor and the donor. Under 
those conditions phase loss can lead to sequential hopping. The interplay between the 
Sequential hopping and tunneling has been adressed by Medvedev and Stuchebrukhov 
|21| . They obtained the expression for rate by forming a dynamic correlation function 
of the couplings or the hopping parameter. Alternate methods to arrive at the rate 
expression include, density matrix formulation of the problem and equation of motion 
for the reduced quantum system interacting [22l [23l l24l [251 l26l [27] . 

Our present work focuses on obtaining an explict expression for transfer rate in 
a system where the redox is solvated and is connected by a molecular chain to a 
electrode. More specifically our attention is concentrated towards obtaining the voltage 
dependecy of the electron transfer rate. As with the previous authors who tackled the 
problem within the time-dependent fluctuation framework, we assume that at t=0 
the electron is in the donor and express the rate of the electron arrival at the final 
continuum states in the electrode using t-matrix. In this way our treatment of the 
problem is different from the previous treatments by other authors, wherein the rate 
expression was derived by expressing the rate as a time dependent correlation function 
and by partial tracing of the density matrix. It should also be noted that all the earlier 
works for obtained transfer rate were mainly oriented towards a DBA system or some 
likewise where both the donor and acceptor were either solvated or continuum of 
states, while in our system considered only the donor is solvated while the acceptor is 
an electrode. 



2 Model and Calculation 



The following Hamiltonian is considered as a model for a system of N-sites (single 
energy level £j) connecting a reservoir or electrode with energy levels tk to a redox 
couple with energy level e r . The electrode and the redox is connected to the N-site 
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chain through the 1 and N- site. The Hamiltonian is a tight binding model hamiltonian 
with additional interactions to account for the interaction of the adsorbate and the 
redox with the polrization modes. This is acheieved by treating these polarization 
modes as classical phonon modes. For ease of clarity the Hamiltonian can be split into 
three parts H e i, H chain and H ph , where the H e i refers to the electrode part of the 
Hamiltonian , H c h a in represents the tight-binding Hamiltonian of the chain + redox 
couple and the H p h is the used for modelling the classical phonon modes. 



U e i = Y, k e k n k + T, k [v kl c[ci + h.c] 



(1) 



H 



ph 



T, v -hw v q v 



(2) 



H. 



chain 



^ ejUj + y^{vjc\ c i+ i + h.c} + (e r 



-1,N 



^ hw u g v q v )n r + {v nr c\ l c r + h.c} 

(3) 



It is clear that the matrix form for the H c h ain is a N+l x N+l tridiagonal matrix 
with only non-zero entries along the main diagonal and sub-diagonals. For complete- 
ness we also give below the matrix form of the hamiltonian 



H 



chain 



o 



/ ei v . 

V €2 V 
V 

Q ... v e N v rn 

H = H e ( + Hchain + Hph 



(4) 



(5) 



To get the current the formalism of Ratner et al is employed. Thus the proce- 
dure is to get the rate for the transition from redox (r) to the electrode (k) under the 
influence of voltage W and then perform a thermal averaging of the transition rate. 

rate = — S fe [l - /(e - eW)] \ T kr | 2 5(e r + Y< u hg v q u - e k ) (6) 



Tfe r = VklGlnVn 



(7) 



pit) 



(- v y 



det[H, 



chain] In 



(8) 
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The Hchain refers to the Hamiltonian of the chain along with the self-energy correc- 
tion arising from the coupling of it's 1-st site with the electrode. In effect, this amounts 
to adding the self-energy to the 1,1 element of the K c hain- The Gi„ is calculated based 
on the following crucial observation. Let d n denote the determinant of n x n matrix 
consisting of e in the diagonal elements and v in the upper and lower sub-diagonal and 
the rest of the elements of the matrix being zero. 

( t. ) " ( i o" 2 ) ( tl ) < 9 » 



By repeated application of the above recursion relation, and noting that do = 1 
and di = e one arrives at the following result 

\rt+l _ \n+l 

d n = \ ? (10) 



a 1i2 = i^EEiZ 



Employing this result in the present Hamiltonian, yields the below given result for 
the determinant of Yl c hain 

det[H cfco . n ]i n = d n +iAd n _i (12) 



Where in the above expression d n is the same as the one defined above . Care 
should be taken to note that the e used in the definition of Ai,2 is to be replaced with 
e — e». This is done since the quantity of interest is the inverse of el — H c hain 



In the above expression the green's function element obtained was of the isolated 
chain Hamiltonian, but the green's function element to be employed in the formal- 
ism should correspond to the full total Hamiltonian. This can be obtained from the 
isolated chain Hamiltonian's green's function by resorting to a Dyson equation with 
the interaction of the electrode and redox treated as a perturbation. The resulting 
expression is as follows: 

GlN = \1 - G° Gl N v 2 G° 1 (13) 



Where in the above expression the isolated green's function elements are given by 
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Gnn ~ d n + t Ad n _ 1 (14) 
^ rr e - e r - Y, v hjJ v g v q v ± iS 



Alternatively one can arrive at the expression for the Gi„ element by considering 
only the Hamiltonian for the chain, }i c h a in and replacing the effect of the reservoir 
and redox by the corresponding self-energies. Working within the Wide-Band approxi- 
mation where the self-energy of the electrode is given by ( — iA ). Both the procedures 
leads to the expression for Gin as 



Where in the above expression wide-band approximation has been employed for 
the self-energy of the reservoir and the self-energy of the redox is assumed to first 
approximation as jr^^S^—^ 



The above expression for Gi„ is to be substituted in the rate expression. The delta 
function can also be split and hence the rate expression takes the form as shown below 

2-7T f 

n{e,q v ) = — / deE fe [l-/(e-eW)] | v kl | 2 S(e-e k ) \ v nr | 2 | Gi r | 2 5(e r +E v hg„q v -e) 

(17) 



Defining ty \ vti \ 2 Sfc(5(e — e^) = A the above expression reduces to the form below 



«(e, q v ) - |A[1 - f(e - eW)] \ G ln | 2 | v nr | 2 (18) 



The £ in the above expression refers to the energy at which the electron transfer 
takes place and hence to get the net current flow it is required to sum all the contri- 
butions from different energy windows. This summing of different contributions can 
be performed by integrating with respect to e. That is a sum over e is replaced with 
J dep(ef) where p(e/) is the density of states at the fermi-surface of the electrode . 
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«(«„) = \±j dep(e f )[l - f(e - eW)] | G ln | 2 | v nr | 2 (19) 
n{q v ) - jUp( C/ ) 1 de[l - /(e - eW)] \ G ln | 2 | v nr | 2 (20) 



It is now required to do thermal averaging which is to integrate the above expres- 
sion for current over all the polarization modes ,q v , with a weighing factor of exp~ l3E . 
Alternatively the same process can be viewed as summing over all possible initial 
states of the redox with respect to polarization modes with a suitable weighing factor 
depending on the energy of the polarization modes. 



The above expression can be evaluated by employing the single reaction co-ordinated 
Q for the polarization modes q v . Resorting to single which in effect can be represented 
as a change of co-ordinates, the crucial observation is that the introduction of single 
reaction co-ordinate amounts to replacing T, l ,huj^g v q l , by 2\q and the quadratic term 
in the exponential for the thermal average is to be replaced by \q 2 + 2Xq. Thus the 
expression for net current takes the form as shown below 

1 2 



A/>(e,)tW | 2 J de[l - f(e - eW)} J d qe -^ 2+2X ^ \ G ln (q) | 2 | (21) 

{-v) 71 - 1 

Gln{q) = ( rj I JT d .) K- 1+ »A rf „- 2 M/ ( 22 ) 



dqe -0(\ q 2 +2\q) (23) 



While the partition function Z is an elementary gaussian integral. The innermost 
dq integral can be evaluated by using partial fractions. For brevity defining D n = 
d n + iAd n _i, the Gi„ takes a form as below. 

{-vY-\e-e r + 2\q) 
ln D n (e~e r + 2\ q )-v^D n ^ [ ' 



The innermost integral over dq can be explictly evaluated 

dqe -P(X q * + 2X q) | |2 = l^L_ [Ae -^(Qo+2A 2 )/4A (e _ ^ + Qo) 2 

+ Be -^(Qo+2A 2 )/4A (e _ £r + q o)2] (25) 
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Where A = e't r (D„_ 1 D„"-D„_ 1 D„) aIld B = 0it r (D„_ 1 D^-D„_ 1 D„) and Q° arG COm - 

plex conjugates of each other. Physically this corresponds to the pole in Gi n (q). 

Qo = v 2 nr ^-(e-e r ) (26) 



2iA(d n - 2 d n - dn-idn-i) 



2iA(d n - 2 d„ - d n -id n -i) 



P.F 



2 



A I v. 



2n-2 



-/3A 



(29) 

(30) 



The above expression for current can be written in a simplified form by resorting 
to a step function approximation to 1 — f(e — eW) . 



"I 

J €f 



TA/J J J J \ " £ (31) 

_jy 2iA(d n -2d n — d n -\d n -i) 



3 Results and Discussion 



Our main focus of attention in this work is oriented towards explaining the dynamic 
interaction of the redox with the solvent modes and it's subsequent effect on the 
transfer rate. The chief reason being the redox energy levels varies with the q and 
hence even though maxima in the transfer rate is expected when the energy level of 
the redox is same as the energy level of the bridge, but due to the coupling of the redox 
with the solvent modes, the value of q coupled with the resonant energy level maynot 
be a thermally favoured one. Thus it's the competition between the energy difference 
between the redox and the bridge and the probability distribution of the initial state 
of the redox that brings out certain interesting results which we focus on this part. 
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Since there exist a number of parameters as input, unless stated specifically, we 
set the fermi surface to zero ,e/ = 0, A = 0.2 and v nr = Vi, for all through the rest of 
the discussion. Fig U shows the variation of the rate with the applied voltage which 
could be understood easily from a naive argument as follows: As the applied voltage 
is increased, the energy levels available in the electrode for the electrons originating 
from the redox decreases. Thus leading to a rate approaching zero as the voltage is 
increased. The above is a normally expected behaviour while the interplay between the 
polarisation modes and the resonance between the redox and the bridge energy levels 
is clearly visible in Fig 3] Analysis reveals that the mode corresponding to q — — 1 
is the thermally favoured mode, and hence when the energy difference between redox 
and bridge sites equals 2A is the condition for energy matching between the redox and 
bridge. This condition is satisfied for the set of values employed in obtaining Fig [4] 
For the set of values used, the resonance between redox and bridge occurs at 2 eV. 
Thus patterns with local maximas and minimas occur around 2 eV in Fig [H while the 
curve falls smoothly and rapidly around the same value in Fig g] . 

Fig [4] is plotted to show the temperature variation of the rate at zero bias. As 
expected, the graph shows a linear behaviour when plotted as a function of (1/T). 
Taking logarithm of the rate equation reveals only three terms dependent on j3, apart 
from uninteresting constant terms. One of the terms is -0.5 log /3 while the other two 
terms have linear dependency on /3. The contribution from log j3 become dominant 
only in the /3 — > limit. Since this corresponds to the high temperature limit , for all 
practical purposes the log variation of rate has a linear dependency on (1/T), a fact 
borne by Fig [4] 



In order to further investigate the properties of transfer rate at resonance condition, 
we plot in Fig [4] the response of the rate- Voltage for various values of vi. It is 
observed from the plots that for small values of the coupling co-effecients , implying a 
weak bonding between the bridge atoms and between the bridge atoms and the redox, 
drastic variation from the normal behaviour is not seen. But as the value of v\ is varied, 
the observed pattern near the resonance point changes. The changes include a constant 
transfer rate around a small range of voltage near the voltage value corresponding to 
the resonance energy, to the rate dropping to zero in some intermediate voltage regime. 

In Fig 21 we plot the variation of the rate as a function of the difference in energy 
between the bridge and the unsolvated redox energy. In particular we wish to highlight 
the behaviour in the regime where the bridge is situated at a higher energy level than 
the fermi level, we first study the variation for a fixed value of A. The crucial fact 
to observe is that as long as the bridge is located energetically higher than the fermi 
surface, the location of the maximas varies with the actual location of energy level of 
the bridge and not just only upon the energy difference between the redox and the 
bridge as one would normally expect. On examining the plots show that the maxima 



8 



occurs for a fixed value of e r corresponding to a particular value of A. Thus in effect 
while the minima is dependent on the difference between the energy levels, the maxima 
shows a strong dependecy on the solvated energy level of the redox. The effect of A on 
the rate profile is shown in Fig [3] 

4 conclusion 

In this article we have considered a model for electron transfer rate from redox to 
an electrode via a bridge of N- atoms. The redox considered is allowed to interact with 
an electrochemical system whose effect in terms of polarization modes are modelled 
as a bath of classical oscillators with solvent co-ordinates q. We obtained the electron 
transfer rate from redox to the electrode and have performed an averaging over all 
possible initial states of the redox. A final voltage dependent expression for the rate 
is obtained as an integral over all possible energy values available for the incoming 
electron. This integral has been performed numerically. 



We have shown some of the interesting pattern which are otherwise not expected in 
a heuristic analysis. These pattern in the voltage dependency have been attributed to 
the competing process namely,resonance tunneling between the redox and the bridge 
and the solvation of the redox. Additionally we have also shown that by suitable alter- 
ing the voltage range and the coupling co-efficients it's possible to block the electron 
transfer in some intermediate voltage ranges. Moreover, the profile of the transfer rate 
with respect to the energy difference between the bridge and the redox exposes a 
radically different response due to the effect of intervening solvent modes. 
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